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Abstract 

The ground state and the finite temperature phase diagrams with respect to 
magnetic configurations are studied systematically for thin magnetic films in terms 
of a classical Heisenberg model including magnetic dipole-dipole interaction and 
uniaxial anisotropy. Simple relations are derived for the occurrence of the various 
phase boundaries between the different regions of the magnetic orientations. In 
particular, the range of the first and second order reorientation phase transitions 
are determined for bi- and trilayers. 
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1 Introduction 



Recent developments of thin film technologies enable the control of the growth of ultrathin 
magnetic films deposited on nonmagnetic substrates. Due to their challenging application 
as high-storage magnetic recording media, much attention ( |Gradmann 1986] , [AllenspcicE| 



et al. 19901 , Pappas et al. 1990| , |1992|. lAllenspach 1994| , |Li et al. 1994i [Berger and Hopster 



1996, Garreau et al. 1996 , Farle et al. 1997 , Gubiotti et al. 1997 ) has been devoted to the 



novel properties of these new structures. From a technological point of view, the study 
of the magnetic phase transitions, in particular, of reorientations of the magnetisation is 
playing a major role. As compared to bulk systems, the presence of surfaces and inter- 
faces leads to an enhancement of the magneto-crystalline anisotropy due to spin-orbit 
coupling. The magneto-crystalline anisotropy often prefers a magnetisation perpendicu- 
lar to the surface, while the magnetic dipole-dipole interaction and the entropy at finite 
temperatures favor an in-plane magnetisation. Consequently, as observed in many Fe or 
Co based ultrathin films QGradmann 198(^ , Allenspach et al. 1990 , Pappas et al. 1990| , 
199^ , [Allenspach T99^ , |Li et al. 1994| , [Berger and Hopster 1996t [Garreau et al. 199q ), a 



reorientation from out-of-plane to an in-plane direction of the magnetisation occurs by 
increasing both the thickness of the film or the temperature. Relativistic first principles 
calculations using the local spin density approximation (LSDA) turned out to be suffi- 



ciently accurate to reproduce the critical thickness of the reorientation ([Szunyogh et al 
[19951 , [T997b| , [Zabloudil et al. 1998[) . 

In the case of thin Ni films deposited on a Cu(OOl) surface, an opposite behavior was 
revealed: the magnetisation was found to be in-plane for less than 7 Ni monolayers, how- 
ever, it became perpendicular to the surface beyond this thickness. Below the switching 
thickness near K, even an increase of the temperature induces a similar reversed reorien- 
tation ( [Farle et al. 1997[ , [Gubiotti et al. 1997| ). The main origin of the above reorientation 
was attributed to a strain induced anisotropy of the inner layers preferring a perpendicular 
magnetisation ( [Hjortstam et al. 1997[ , [Diberacker et al. 1999[ ). 



Subsequent to the pioneering works of [Mills (1989[) , who predicted the existence of a canted 
non-coUinear ground state for a semi-infinite ferromagnetic system, and that of Pescia and 
Pokrovsky (1990[) who, by using a renormalization group treatment of a continuum vector 
field model, for the first time described the temperature induced (normal) reorientation 
phase transition, a considerable amount of theoretical attempts, mostly by means of 
different statistical spin models, was suggested in order to explain the above findings 
for the thickness and temperature driven reorientation transitions. In some attempts. 
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classical vector spin models were used within the mean field approximation ([Taylor and 
yorfty 1993| , |Hucht and Usadel 1996| , |T99l , [T999a| ,[b|, pO0| , [Jensen and Bennemann T998| , 



Hu et al. 1999[) or in terms of Monte Carlo simulations ([Taylor and Gyorft'y 1993] , ^erena 
St ai 1993[ , [Chui 1995| , [Hucht et a/.T995t [Hucht and Usadel 1996| , [Maclsaac et al. 1996[) . 
A quantum-spin description of reorientation transitions has also been provided in terms 
of spin-wave theory ( Pruno 1991[ ), mean field theory ([Moschel and Usadel 199"4[ , [1995[) , 
many-body Green function techniques ( [Frobrich et al. 2000a| ,[E[, [Jensen et al. 2000[ ), and by 
using Schwinger bosonization ( [Timm and Jensen 2000[) . Although, the mean field theory 
is not expected to give a sufficiently accurate description of low-dimensional systems, it 
turned out, that it is a successful tool to study spin reorientation transitions and yields 
qualitatively correct predictions ([Moschel and Usadel 1994[ , [1995[ , [Hucht and Usadel 1997[ , 
1999a POOOI , [Jensen and Bennemann 199S , Hu et al. 1999[) . It also should be noted that 



an itinerant electron Hubbard model revealed the sensitivity of reorientation transitions 
with respect to electron correlation effects ( [Herrmann et al. 1998[ ). 

For layered systems, the following simple model Hamiltonian can be used to study reori- 
entation transitions, see e.g. ( [Taylor and Gyorff'y 1993[) , 
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where Spi (| Spj |= 1) is a classical vector spin at lattice position i in layer p and Vpi^gj is a 
vector pointing from site {p, i) to site {q,j) measured in units of the two-dimensional (2D) 
lattice constant of the system, a. Although our previous calculations of the Heisenberg ex- 
change parameters in thin Fe, Co and Ni films on Cu(OOl) showed some layer-dependence 
( |Szunyogh and Udvardi 1998[ , [1999[ ) , in the first term of Eq. (|I|) we only consider a uniform 
nearest neighbor (NN) coupling parameter J throughout the film. Similarly, as we ne- 
glect the well-known surface/interface induced enhancement of the spin-moments, we use 
a single parameter u = HqIj? /Ana^ (with /iq the magnetic permeability and n an average 
magnitude of the spin-moments), characterizing the magnetic dipole-dipole interaction 
strength in the third term of Eq. (|l]). As revealed also by first principles calculations - 
see [Weinberger and Szunyogh (200C ) and references therein - , the uniaxial magneto- 
crystalline anisotropy depends very sensitively on the type of the surface/interface, the 
layer-wise resolution of which can vary from system to system. Therefore, the correspond- 
ing parameters, Ap, in the second term of Eq. (0) remain layer dependent: the variety of 
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these anisotropy parameters leads to rich magnetic phase diagrams covering the experi- 
mentally detected features mentioned above. For example, in a previous study QUdvardi 
et al. 1998D we pointed out that, even in the absence of a fourth order anisotropy term, for 
a very asymmetric distribution of Ap with respect to the layers, the Heisenberg model in 
Eq. (|I]) can yield a canted (non-collinear) ground state. This feature cannot be obtained 
within a phenomenological single domain picture. 

As what follows, we first investigate the possible ground states of the Hamiltonian, E q. 
Then we perform a systematic mean field study of the different kind of temperature in- 
duced reorientation transitions, devoting special attention to the case of bi- and trilayers. 
Specifically, we define general conditions for the reversed reorientation. Most authors in 
the past focused on proving the existence of different reorientations and detected only 
some parts of the phase diagram, where first and second order phase transitions occurred. 
In here, we describe the full range of uniaxial anisotropics (Ap), for which first or second 
order reorientation phase transitions can exist. Finally, we attempt to summarize the 
results and impacts of a mean field approach. 



2 Ground state 

Confining ourselves to spin-states in which the spins are parallel in a given layer, but 
their orientations may differ from layer to layer, i.e. 



Spi = Sp = (sin(6'p) cos(0p), sin(6'p) sin(0p), cos(6'p)) 



(2) 



where Op and 0p are the usual azimuthal and polar angles with the z axis normal to the 
planes, the energy of layers per 2D unit cell can be written as 
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with Upq being the number of nearest neighbors in layer g of a site in layer p, and Apq the 



magnetic dipole-dipole coupling constants, see Appendix in QSzunyogh et al. 1995|) 
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which is vahd for square and hexagonal 2D lattices, such as the (100) and (111) surfaces 
of cubic systems. (/ the three dimensional unit matrix, while ® stands for the tensorial 
product of two vectors.) For three-dimensional (3D) translational invariant underlying 
parent lattices ( Weinberger 1997 ), Apg depends only on \p ~ q\, i.e. the distance between 
layers p and q. In table ^we summarize these constants for the first few layers of the most 
common cubic structures. As the magnetic dipole-dipole interaction is clearly dominated 
by the positive first (and second) layer couplings, for ferromagnetic systems (J > 0), a 
minimum of the energy in Eq. (^ corresponds to the case when all polar angles, 0p, are 
identical. Therefore, the 0-dependence in Eq. (H) disappears and the expression, 
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has to be minimized with respect to the 6p. The corresponding Euler-Lagrange equations 
are then 
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Obviously, a uniform in-plane {{6p = vr/2}) and a normal-to-plane {{6p = 0}) orientation 
satisfy Eq. (|^). The energies of these two particular spin-states coincide, if 
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which defines an (A^—1) -dimensional hyper-plane in the A^-dimensional space of param- 
eters {Xp/uo}. If the magnetisation changes continuously across this plane, in its vicinity 
there should exist solutions with canted magnetisation. Moreover, the saddle points of 
the energy functional in Eq. (^, 
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define the boundaries of the canted zone, 
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for the uniform in-plane and normal-to-plane magnetisations, respectively. For a bilayer 
we derived explicit expressions of Eqs. (H) and (p!OD, see [Udvardi et al. (1998|) . 



In order to study the canted region, instead of solving the Euler-Lagrange equations. 



Eq. (|^), directly, we fixed a configuration O^.O^ 
demanding that Eq. (^) must be satisfied, i.e. 



6*^ and determined Ai,A2...AAr by 
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Substituting these parameters into Eq. (||), one easily can express the difference of the 
energies between the corresponding configurations as. 
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Note that the position of the minimum {6*} as well as the minimum energy E]\f{{6*}) 
are functions of the parameters J, uj, and {Ap}. Obviously, whenever the parameters fall 
into the region between the two hyper-planes defined by Eqs. (|^) and (^^), the energy 
of non-collinear states is always below or equal to the energy of the collinear in-plane or 
normal-to-plane solutions. 

We showed that for a bilayer, Ai = A2 = 3(y4ii + Ai2)uj /A implies a collinear ground state 
spin configuration ( [Udvardi et al. 1998|) . This state is, however, continuously degenerate, 
i.e. the energy is independent of the orientation of the magnetisation. Such a critical 
point in the phase diagram also exists for multilayers (A^ > 3). Namely, from Eqs. (|12D 
and (|T3D it follows that for = ■ ■ = = 0* , ^iv({^;}) is independent of e\ In 
terms of Eq. (0), the corresponding point in the parameter space {Xp/uj} is given by 
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Evidently, the hyper-planes given by Eq. and Eq. (|10|) touch the hyper-plane, separat- 
ing the in-plane and normal-to-plane regions, Eq. (^, at the point defined by Eq. (|14D- 
It is worthwhile to mention that this is the only point where canted collinear solutions 
can exist. This critical point was also found by Hucht and Usadel (199(: ) for a monolayer, 
but they did not prove its existence for multilayers. 



3 Finite temperature 

Introducing the following coupling constants 

Cpg = npqJ + ^ujApg and c^^ = UpgJ - uApq , (15) 
the molecular field corresponding to the Hamiltonian Eq. (0) at layer p can be written as 

N 

HP {6, (P) = -J2 Cpg [nfp sin {9) cos (0) + sin {9) sin (0)] 

(IS) 

-Zl^p^K^os (^) - Apcos^ {9) , 

9=1 

where = (s") (a = x, z). Similar to the ground state (see Section II.), due to the in- 
plane rotational symmetry of the above effective Hamiltonian, the in-plane projections of 
all the average magnetic moments nip are aligned. Therefore, by choosing an appropriate 
coordinate system, can be taken to be zero in Eq. ([I6|). The partition function is then 
defined by 

N 

Z = l[Zp, (17) 
p=i 



Zp = 2Ti j exp{/3[6^cos(0) + ApCos^(^)]} Jo(-2/56pSin(^))sin(^)c^^, (18) 

7r 

2 

where 

N 

bf^ = T.4'W'^^ (19) 

9=1 

f3 = l/lksT), ks the Boltzmann-constant and T the temperature. The minimization of 
the free-energy with respect to the average magnetisations leads to the following set of 
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nonlinear equations 
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In Eqs. (|T^, (13) ^^'i (|2l|) , Jq and Ji denote Bessel functions of zero and first order, 
respectively ( [Abramowitz and Stegun 1972|) . 

By using a high temperatures expansion, Eqs. (|20|) and (^iD become decoupled (see Ap- 
pendix). Consequently, the magnetisation can go to zero either via an in-plane or via a 
normal-to-plane direction at temperatures and Tj, respectively, and the higher one of 
them can be associated with the Curie temperature Tc- Clearly, an out-of-plane to in- 
plane reorientation phase transition can occur only when the ground state magnetisation 
is out-of-plane and < = Tc- In the case of a reversed reorientation transition, the 
ground state magnetisation has to be in-plane (or canted) and < = Tc- 

Expanding T^ and up to first order of the anisotropy parameters Xp, leads to the 
following expressions (see Appendix) 
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The above expressions imply that, if the anisotropy parameters Xp are small, T^ is larger 
than Tz- As the anisotropy parameters are increasing, the difference between and T^ 
decreases. The two temperatures coincide, if the following condition is fulfilled. 
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Above the hyper-plane determined by Eq. (|2J), i.e. for T^. < T^, the uniaxial anisotropy is 
large enough to keep the magnetisation normal to the surface as long as the temperature 
reaches Tq- 

First principles calculations on (Fe,Co,Ni)/Cu(001) overlayers revealed ( |Ujfalussy et at 



L99(j| , [Szunyogh et at. 1997a , Szunyogh and Udvardi 1998 , 1999| , Uiberacker et at. 1999 ), 



that the uniaxial magnetic anisotropy energy and the magnetic dipole-dipole interaction 
are two to three orders of magnitude smaller than the exchange coupling. Thus, for 
physically relevant parameters, the boundaries of the canted ground state fixed by Eqs. 
and (^) are close to the hyper-plane defined by Eq. (0). Apart form this tiny range of 
canted ground states, temperature induced out-of-plane to in-plane reorientation can 
occur in the parameter space {Xp/u} between the two hyper-planes given by Eqs. and 
p^). It is worthwhile to mention that the positions of these hyper-planes are determined 
only by the magnetic dipole-dipole constants Apg. 

An example for an out-of-plane to in-plane reorientation transition in a 5-layer thick 
film is shown in figure |l]. Neglecting the fourth order anisotropy terms, the parameters of 
the system have been chosen identical to those characteristic to a Co5/Au(lll) overlayer 
( [Udvardi et at. 1998[ ). Due to the highly asymmetric distribution of the Xp with respect 
to the layers, the system has a non-collinear canted ground state. As the temperature 
increases, the magnetisation in each layer turns into the plane of the film. The system 
keeps its non-collinear configuration up to the reorientation transition temperature (~ 
0.9 J/ ks), above which it is uniformly magnetized in-plane up to the Curie temperature 
(~ 3.8J/kB). 

The temperature induced reversed reorientation transition, found experimentally in Ni„ / 
Cu(OOl) films for n < 7, has successfully been described by Pucht and Usadel (1997|) , who 
used a perturbative mean field approach to the model given in Eq. (|l|). Using the same 
parameters, we solved the mean field equations ( pO] ) and ( pT] ) and reproduced the reversed 
reorientation transition without any perturbative treatment. The results for a 4-layer film 
are shown in figures ^ and |^. Although the distribution of the anisotropy parameters is 
asymmetric, the calculation resulted in identical magnetisations in the first and fourth 
layer as well as in the second and third layer. Moreover, the angles of the magnetisations 
in the different layers are almost identical: in the whole temperature range the largest 
deviation is smaller than 6 x 10~^ rad. 
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For a bilayer {N = 2) the hyper-planes (0) and (^41) reduce to the hnes 



^ + ^ = Ua, + A,,) and ^ + ^ = ^An + ^,), (25) 
io UJ 2 UJ UJ 2 

respectively. Apparently, the two lines do not intersect. As a consequence, a reversed 

reorientation can occur only if the number of layers in the film exceeds two. The same 

conclusion has been drawn by [Hucht and Usadel (1997| ) using a perturbative treatment of 

the anisotropy parameters. Nevertheless, it is interesting to note, that the region in the 

parameter space {Xp/uj} of canted ground states, bounded by the lines defined by Eqs. (H) 

and (0), always overlaps the region, where the magnetisation goes to zero via in-plane 

orientation. Thus, in this overlapping region an out-of-plane (canted) to in-plane, i.e. 

reversed reorientation transition can indeed occur. The corresponding parameters, u and 

{Ap}, are, however, most likely beyond the physically relevant regime. 

For a bilayer, in figure § the different regions of phase transitions in the respective pa- 
rameter space are shown. In regions I and V there is no temperature driven reorientation 
transition and the magnetisation remains in-plane and normal-to-plane, respectively, un- 
til Tc is reached. In the regions II and III, the magnetisation turns into the plane from 
a canted or a normal-to-plane ground state, respectively. As discussed above, in region 
IV, a reversed reorientation can occur from a canted ground state to a normal-to-plane 
direction. 

The order of the reorientation transition at finite temperatures has been studied in the 
literature by mean field and Monte Carlo methods. Most authors concluded ( [Hucht et al 



[I995| , Pucht and Usadel 1996| , [Maclsaac et al. 1996| ) that the reorientation transition in 



a monolayer is of first order. For a bilayer, within the mean field approach, a relatively 
small range in the vicinity of Ai = A2 was found, where the system underwent a first 
order reorientation transition ([Hucht and Usadel 1996|) . In what follows, we establish a 



simple, general criterion for the order of the reorientation transition. Suppose that the 
ground state magnetisation is in-plane and its normal-to-plane component appears at 
the temperature Trz- Since near T^z the ^-component of the magnetisation is small, the 
exponential function in Eq. (^) can be expanded up to first order in m^, leading to the 
homogeneous linear equations 

N 

9=1 

where 

Cpq — ^pq ~ PrzCpqZp , (27) 
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O-rr f 

Zp = — cos\6) exp [przXp cos{9f] M-iPrzbp sm{9)) sm{9)de , (28) 
Zp J 

and Prz = ^/{kBTrz)- Note, that the Cp^ depend only on the in-plane component of the 
magnetisations m^', which has to satisfy Eq. (^) for the case of = 0. Eq. (pB]) has a 
non-trivial solution only, if the determinant of the matrix = {C'pg} is zero. Evidently, 
this is the condition which determines T^z- Similarly, one can easily find the corresponding 
equation for T^x, where the in-plane component of the magnetisation appears in a normal- 
to-plane spin configuration, 

N 

E^>^ = ' (29) 

9=1 

with 

(31) 

and Prx = ^/iksTrx)- It is easy to show, that Eqs. ( p6D and ( [29| ) directly follow from a 
stability analysis of the mean field free energy in the vicinity, where the corresponding 
components of the magnetisations vanish. 

The mean field equations, Eqs. (pOD and (|2T|), always have an in-plane and a normal-to- 
plane solution with magnetisations ml ^ 0, ml = and ml ^ 0, = 0, respectively. 
Between Trx and T^z, a canted solution can exist with m^ ^ and ml ^ 0. Among the 
above three phases, the physical one belongs to that, which has the lowest free-energy. 
In figure §a, the free-energy of a system possessing a second order normal-to-plane to 
in-plane reorientation transition is schematically shown. The ground state magnetisation 
is perpendicular to the surface of the substrate. At T^x an in-plane component appears in 
the magnetisation. The normal-to-plane component of the magnetisation vanishes at the 
temperature Trz{> Trx)- A similar picture for a first order transition is shown in figure 
Obviously, one can conclude that, if T^-x < Trz, a second order normal-to-plane to 
in-plane reorientation phase transition occurs, whereas, if T^z < Trx, the reorientation 
transition is of first order. In the case of a reversed reorientation, the relation between 
Trx and Trz is just the opposite as before: a second order transition occurs, if Trz < Trx, 
while for Trx < Trz the transition is of first order. At the boundary of the regions, where 
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second order and first order phase transitions occur, the two temperatures, T^x and T^z-, 
must evidently coincide. 

In figure ^, the region of first order reorientations (III F) and that of second order reori- 
entations (III S) are shown in the phase diagram for a bilayer. Note, that figure |^ in fact 
represents figure ^ on an enlarged scale for the parameters < \i^2/^ < 18. This picture 
is consistent with the observation of Hucht and Usadel (1996 ) for the range of the first 
order reorientation phase transitions, as they performed investigations very close to the 
critical point only. In that case, by keeping Ai + A2 fixed, figure ^ implies a very narrow 
range for the first order transitions. 

The phase diagram of the trilayer case (A^ = 3) is shown in figure |^. Apparently, the same 
regimes exist as in the case of a bilayer. The region of first order reorientation transition 
forms now a 'sack', touching the plane defined by Eq. (^ at the critical point given by 
Eq. ([T^). The sack is covered by the plane separating the area where normal-to-plane 
to in-plane reorientation occur and the area, where the magnetisation remains normal- 
to-plane up to the Curie temperature, see Eq. (0). The regime of reversed reorientation 
transitions, part of the region of second order transitions, is, however, out of the segment 
of the parameter space depicted in figure 0. 

Numerical calculations using different magnetic dipole-dipole coupling strengths u (see, 
in particular, figure ^ yields practically the same boundaries in the {Xp/u} parameter 
space of the phase diagrams for both the bilayer and trilayer cases. The only exception 
is the region of the canted ground states, which rapidly opens up with increasing u. The 
established universality of the phase boundaries nicely confirms that the reorientation 
phase transitions, as long as u gets comparable to J, are a consequence of the competition 
between the uniaxial anisotropy and the magnetic dipole-dipole interaction. 



4 Conclusions 

In the present paper we provided a full account of the ground states and of the finite 
temperature behavior of a ferromagnetic film of finite number of layers, as described by 
the classical vector spin Hamiltonian, Eq. (|l|), including exchange coupling interaction, 
uniaxial magneto-crystalline anisotropics and magnetic dipole-dipole interaction. We 
derived explicit expressions for the boundaries of the regions related to normal-to-plane, 
canted and in-plane ground states in the corresponding parameter space. We concluded 
that within the model, defined by Eq. (||), canted ground states are ultimately connected 
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to non-collinear spin-configurations. In addition - so far establislied for monolayers only 
QHuclit and Usadel 1996| ) - for any thickness of the film we proved the existence of a 
critical point, where the ground state energy of the system is independent from a uniform 
orientation of the magnetisation. 

We also investigated intensively the finite temperature behavior of the system in terms of 
a mean field theory. By using a high temperature expansion technique, we showed that 
the Curie temperature of a ferromagnetic film can be calculated by solving an eigenvalue 
problem, which, for the case of a bulk system and by neglecting anisotropy effects, leads 
to the well-known expression of Tc- The main part of the present study has been de- 
voted to the reorientation phase transitions, which play a central role for applications of 
thin film and multilayer systems as high-storage magnetic recording devices. Both the 
normal-to-plane to in-plane and the in-plane to normal-to-plane (reversed) temperature 
induced reorientation transitions have been discussed and the corresponding regions in 
the parameter space have been explicitly determined. In accordance to previous studies 
QHucht and Usadel 1997|) , we showed that, for physically relevant parameters, reversed 
reorientation can occur only for films containing three or more atomic layers. By investi- 
gating the order of reorientation phase transitions, we found well-defined conditions for 
the first and the second order phase transitions and presented the corresponding regions 
for bi- and trilayers in the respective parameter spaces. 

In conclusion, we have shown that a mean field treatment of a classical vector spin model 
recovers most of the important phenomena observed in magnetic thin film measurements 
at finite temperatures. Without any doubt, due to the lack of mean field theories for low- 
dimensional systems, some of them have to be refined by using more sophisticated methods 
of statistical physics (see Introduction). In particular, for very thin films (monolayers), 
the mean field theory predicts a Tc much higher than the random phase approximation 
(RPA). However, by rescaling the temperature, the orientations of the magnetisation 
become fairly similar in both approaches ( [t'robrich et al. 2000a| ,[5D. As far as the first 
principles attempts ( |Szunyogh et al. 1995| , |1997b| , [Szunyogh and Udvardi 1998| , |1999| , 
[Uiberacker et al. 1999| , Pa] da et al. 2000D are concerned, which are currently able to 
calculate realistic parameters for a model hke Eq. (|ID, the technique, presented and applied 
here, provides a simple and quick tool to study the finite temperature behavior of thin 
magnetic films. As the measurements are performed at finite temperatures, while first 
principles calculations usually refer to the ground state, such a procedure would improve 
the predictive power of ab-initio theories. It also should be mentioned, that first attempts 
to an ab-initio type description of thin magnetic films at finite temperatures, i.e. taking 
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into account the coupling of the itinerant nature of the electrons and the spin degree of 
freedom, are currently under progress ( [Razee et al .1) . 
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Appendix: Derivation of the Curie temperature 

In the high temperature limit {(3 — >■ 0, jdhp^^"^ — > 0) the partition function as given by Eq. 
(|T^) can be written up to the first order of the magnetisation as 



1 



1 




(32) 



-1 



-1 



Similarly, for the magnetisation in Eq. (^I]) the following approach can be used 



1 




-1 




(33) 



where an external magnetic field H has been added to the Hamiltonian in Eq. ([TB|). By 
substituting the expansion 




into Eq. (^3]) follows 




(35) 
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Requiring non-vanishing magnetisation at zero external field results in the following eigen- 
value problem 



N 
q=l 



'pq 
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m„ 



(36) 



Let Tz denote the highest value of T, for which Eq. ( p6D is satisfied, i.e. above which no 
spontaneous normal-to-plane magnetisation can exist. A similar procedure can be applied 
in order to determine T^, that is, the temperature, at which the in-plane magnetisation 
vanishes. Quite obviously, by neglecting anisotropy effects, for a bulk system the Curie 



temperature Tc = 



Tz is given by the well-known formula 
nJ 



T, 



c 



3kf 



(37) 



where n denotes the number of nearest neighbors in the bulk. 

With exception of very open surfaces such as the BCC(lll) surface (see table |lD, the 
magnetic dipole-dipole coupling constants Apg fall off exponentially with the distance 
between layer p and layer q. Therefore, as an approximation we neglect all Apg for |p-q'|> 1, 
which, by recalling the nearest neighbor approximation for the exchange coupling, implies 
that the matrix formed by the elements c^^ is tridiagonal. The non-vanishing elements 
are then written as 



"pp 



nil J - Aiiu and Cpp_i = Cp_ip = n^J - Ayi'^ . 
Setting Ap = 0, the solution of the eigenvalue problem Eq. ( pSj ) yields 

nil J — ^111-^ + 2(ni2 J — ^121-1^) cos 



2^(0) ^ ±_ 



N- 



(St 



(39) 



with the components of the corresponding normalized eigenvector Up = ^^2/ {N + 1) 
X sm {P7t/{N + 1)). Substituting T^^ into Eq. (||) and using first order perturbation 
theory with respect to Xp, one gets Eq. (p3D for T^. Again, a similar procedure applies for 
deriving in Eq. (53). 
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structure 


An 


Al2 


Ai3 


Au 


SC (100) 


9.0336 


-0.3275 


-0.00055 


< 10-5 


SC (111) 


11.0342 


5.9676 


0.4056 


0.0146 


BCC (100) 


9.0336 


4.1764 


-0.32746 


0.01238 


BCC (111) 


11.0342 


15.8147 


5.9676 


-4.0662 


FCC (100) 


9.0336 


1.4294 


-0.0226 


0.00026 


FCC (111) 


11.0342 


0.4056 


0.00113 


< 10-5 



Table 1: Dipole-dipole coupling constants as defined in Eq. (|) for surfaces of cubic 
structures. 
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Figure 1: Out-of-plane to in-plane reorientation transition in a 5- layer system (Ai/J = 
0.26, A2,...,5 — 0, uj/ J — 0.0056). The z- and the x-components of the average magnetisa- 
tion as well as the angles of the magnetisation with respect to the normal of the surface 
are shown in the upper and the lower panels, respectively. 
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Figure 2: Normal-to-plane and in-plane components of the layer resolved magnetisation 
for a film of 4 atomic layers exhibiting a reversed reorientation transition. The parameters, 
representative to Ni were taken from [Hucht and Usadel (1997|) : J = 1, Ai/ J = — 3.5x 10~^, 
\i/J = 1.5x10"^ {i > 1), uj/J = 5xl0~^. The inset shows the vicinity of the reorientation 
transition on an enlarged scale. 
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Figure 3: Variation of the angle of the average magnetisation for the system in Fig. ^ 
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Figure 4: Phase diagram of the magnetic ground states and the reorientation phase 
transitions for a bilayer. The magnetic dipole-dipole couphng strength cu — O.OIJ. I: 
in-plane magnetisation up to Tq; II: canted ground state with reorientation transition to 
in-plane orientation; III: normal-to-plane ground state with reorientation transition to 
in-plane orientation; IV: canted ground state with reversed reorientation transition; V: 
normal-to-plane magnetisation up to Tq- 
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Temperature Temperature 



Figure 5: Schematic picture of the free-energy in the case of a second order (a) and a 
first order (6) reorientation phase transition. As indicated by arrows, the sohd, dashed 
and dotted hues refer to the in-plane, normal-to-plane and canted mean field solutions, 
respectively. 




Figure 6: Phase diagram of reorientation transitions for a bilayer (A^ = 2). Labels I, II, 
III and V denote the same regions as in figure ^; regime III, however, is partitioned into 
regions referring to first (III F) and second (III S) order normal reorientation phase tran- 
sitions. The corresponding boundary lines are shown for different values of the magnetic 
dipole-dipole interaction strength uo, measured in units of J. 
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Figure 7: Mean-field phase diagram of reorientation transitions for a trilayer (A^ = 3). 
I: in-plane magnetisation up to Tc, III S: second order normal-to-plane to in-plane 
reorientation, III F: first order normal-to-plane to in-plane reorientation, V: normal-to- 
plane magnetisation up to Tc- 
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